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ABSTRACT 

We examine the orbits of satellite galaxies identified in a suite of N-body/gasdynamical sim- 
ulations of the formation of L* galaxies in a ACDM universe. The numerical resolution of 
the simulations allows us to track in detail the orbits of the ~ ten brightest satellites around 
each primary. Most satellites follow conventional orbits; after turning around, they accrete 
into their host halo and settle on orbits whose apocentric radii are steadily eroded by dynam- 
ical friction. As a result, satellites associated with the primary are typically found within its 
virial radius, r v ; r , and have velocities consistent with a Gaussian distribution with mild radial 
anisotropy. However, a number of outliers are also present. We find that a surprising number 
(about one-third) of satellites identified at z = are on unorthodox orbits, with apocenters 
that exceed their turnaround radii. These include a number of objects with extreme velocities 
and apocentric radii at times exceeding ^3.5 r v i r (or, e.g., > 1 Mpc when scaled to the Milky 
Way). This population of satellites on extreme orbits consists typically of the faint member 
of a satellite pair whose kinship is severed by the tidal field of the primary during first ap- 
proach. Under the right circumstances, the heavier member of the pair remains bound to the 
primary, whilst the lighter companion is ejected onto a highly-energetic orbit. Since the con- 
current accretion of multiple satellite systems is a defining feature of hierarchical models of 
galaxy formation, a fairly robust prediction of this scenario is that at least some of these ex- 
treme objects should be present in the Local Group. We speculate that this three-body ejection 
mechanism may be the origin of (i) some of the newly discovered high-speed satellites around 
M31 (such as Andromeda XIV); (ii) some of the distant fast-receding Local Group members, 
such as Leo I; and (iii) the oddly isolated dwarf spheroidals Cetus and Tucana in the outskirts 
of the Local Group. Our results suggest that care must be exercised when using the orbits of 
the most weakly bound satellites to place constraints on the total mass of the Local Group. 

Key words: galaxies: haloes - galaxies: formation - galaxies: evolution - galaxies: kinematics 
and dynamics. 



1 INTRODUCTION 

The study of Local Group satellite galaxies has been revo- 
lutionized by digital imaging surveys of large areas of the 
sky. More than a dozen new sat ellites have been discov- 
ered in the past c o uple of yea rs i Zucker et al. 20041 
Willman et alTl2005l; Martin et alj 120061: iBelokurov et alj 



2006 



2006, 



20071 : llrwin et alj 120071 ; iMaiewski etail l2007h . due in large 
part to the completion of the Sl oan Digital Sky Survey 
i lYork et alj |2000| : IStrauss et ail l2002h and to concerted cam- 
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paigns designed to image in detail the Andromeda galaxy 
and i ts immediate surrounding s (Ibata et al JfeOOlt iFer guson et al. 
[ 20021; iReitzel & Guhathakurt j |2002|. iMcConnachie et alJ I 2603F 
[ Rich et alj |2004 iGuhathakurta etal] 120061 : iGilbert et alj 120061: 
Chap man et alj|2006l Ibata et al. 2007 submitted). The newly dis- 
covered satellites have extended the faint-end of the galaxy lumi- 
nosity function down to roughly ~ 10 3 Lq, and are likely to pro- 
vide important constraints regarding the mechanisms responsible 
for "lighting up" the baryons in low-mass halos. These, in turn, 
will serve to validate (or falsify) the various theoretical models at- 
tempting to reconcile the wealth of "substructure" predicted in cold 
dark matter (CDM) halos with the scarcity of luminous satellites in 
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the Local Group (see, e.g. Klypin et alj 19991: Bullock et alj 



Benson et a 1. 2002; Stoehr et al 



Kravtsov et alj2004lp"enarrubia et al.ll2007h~ 



2002; Kazantzidis et al. 



2000 



2004 



At the same time, once velocities and distances are secured 
for the newly-discovered satellites, dynamical studies of the total 
mass and spatial extent of the Local Gr oup will gain new impe 



tus. These studies have a long history ([Little & Trem aine 



Zaritsky et al. 1989; Kochanek 1996; Wilkinson & Evans 



1987 



1999 



Evans & Wilkinson! l2000l : iBattaglia et al.ll2005h . but their results 



have traditionally been regarded as tentative rather than conclu- 
sive, particularly because of the small number of objects involved, 
as well as the sensitivity of the results to the inclusion (or omis- 
sion) of one or two objects with large velocities and/or dist ances 
dZaritskv et ai]|l989l ; lKo"chaneklll996l ; ISakamoto et~aT]|2003h . An 
enlarged satellite sample will likely make the conclusions of satel- 
lite dynamical studies more compelling and robust. 

To this end, most theoretical work typically assumes that 
satellites are in equilibrium, and use crafty techniques to over- 
come the limitations of small-N statisti cs when applying Jeans' 
equations to estima t e mas s es (see, e.g., iLittle & Tremainelll987l : 
IWilkinson & Evans! 1 19991 ; lEvans & Wilkinson! l200d) . With in- 
creased sample size, however, follow enhanced opportunities to 
discover satellites on unlikely orbits; i.e., dynamical "outliers" that 
may challenge the expectations of simple-minded models of satel- 
lite formation and evolution. It is important to clarify the origin 
of such systems, given their disproportionate weight in mass esti- 
mates. 

One issue to consider is that the assumption of equilibrium 
must break down when considering outliers in phase space. This 
is because the finite age of the Universe places an upper limit to 
the orbital period of satellites observed in the Local Group; high- 
speed satellites have typically large apocenters and long orbital pe- 
riods, implying that they cannot be dynamically well-mixed and 
casting doubts on the applicability of Jeans' theorem-inspired anal- 
ysis tools. 

To make progress, one possi bility is to explore variants of the 
standard se c ondary infall model dGunn & Got3ll972l; iGotdl 19751 : 

ichf - 



lGunnlll97l [Fmmore & Go ldreich 1984), where satellites are as 



sumed to recede initially with the universal expansion, before turn- 
ing around and collapsing onto the pri mary due to its gravitatio nal 
pull. This is the approach adopted by Zaritsky & White! d 19941) in 
order to interpret statistically the kinematics of observed satellite 
samples without assuming well-mixed orbits and taking into ac- 
count the proper timing and phase of the accretion process. 



sys tenfl r vlr ; the turn around radius is of order r- ta ~ 3 r V ir (see, 
e.g. lWhite etail 1993b . 

We note a few consequences of this model, (a) Satellites out- 
side the virial radius are on their first approach to the system and 
thus have not yet been inside r v i r - (b) Satellites inside the virial ra- 
dius have apocentric radii that typically do not exceed r v i r . (c) The 
farther the turnaround radius the longer it takes for a satellite to turn 
around and accrete and the higher its orbital energy, (d) Satellites 
with extreme velocities will, in general, be those completing their 
first orbit around the primary. Velocities will be maximal near the 
center, where satellites may reach speeds as high as ~ 3 VW- (e) 
Since all satellites associated with the primary are bound (other- 
wise they would not have turned around and collapsed under the 
gravitational pull of the primary), the velocity of the highest-speed 
satellite may be used to estimate a lower limit to the escape veloc- 
ity at its location and, thus, a lower bound to the total mass of the 
system. 

Hierarchical galaxy formation models, such as the current 
ACDM paradigm, suggest further complexity in this picture. 
Firstly, although numerical simulations show that the sequence 
of expansion, turnaround and accretion of satellites described 
above is more or less preserved in hierarchi cal models, the evo- 
lution is far from spherically symmetric (Navarro et al. 
Ghignaetal. 1998; Jing & Sutoll2002l ; rBailin & Steinmetz 



1994 



2005 



Knebe & Wi eBner 200^). Much of the mass (as well as many of 



the satellites) is accreted through filaments of matter embedde d 
within sheets of matter formation (see, e.g., iNavarro et alj|2004h . 
The anisotropic collapse pattern onto a primary implies that the 
turnaround "surface" won't be spherical and that the virial radius 
may not contain all satellites th at have completed at least one or- 
bit ar ound the primary (see, e.g. jBalogh et al.l2000l ; lDiemand et all 
l2007h . 

More importantly for the purposes of this paper, in hierar- 
chical models galaxy systems are assembled by collecting smaller 
systems which themselves, in turn, were assembled out of smaller 
units. This implies that satellites will in general not be accreted in 
isolation, but frequently as part of larger structures containing mul- 
tiple systems. This allows for complex many-body interactions to 
take place during approach to the primary that may result in sub- 
stantial modification to the orbits of accreted satellites. 

We address this issue in this contribution using N- 
body/gasdynamical simulations of galaxy formation in the current 
ACDM paradigm. We introduce briefly the simulations in §[2j and 
analyze and discuss them in § [3] We speculate on possible appli- 
cations to the Local Group satellite population in $4] and conclude 
with a brief summary in §|5] 



In the secondary infall accretion sequence, satellites ini- 
tially farther away accrete later, after turning around from larger 
turnaround radii. The turn-around radius grows with time, at a rate 
the depends on the mass of the primary and its environment, as well 
as on the cosmological model. Three distinct region s surround a 
system formed by spherical s econdary infall (see, e.g.. lBertsc hinger 
1 19851 ; In avarro & Whitdll993b : (i) an outer region beyond the cur- 
rent turnaround radius where satellites are still expanding away, (ii) 
an intermediate region containing satellites that are approaching the 
primary for the first time, and (iii) an inner, "virialized" region con- 
taining all satellites that have turned around at earlier times and are 
still orbiting around the primary. To good approximation, the latter 
region is delineated roughly by the conventional virial radius of a 



2 THE NUMERICAL SIMULATIONS 

We identify satellite galaxies in a suite of eight simulations of 
the formation of L* galaxies in the ACDM scenario. This se- 
ries has been presented by Abadi, Navarro & Steinmetz (2006), 
and follow the same numerical scheme originally introduced by 

1 We define the virial radius, r vlr , of a system as the radius of a sphere of 
mean density ~ A v ; r (z) times the critical density for closure. This defini- 
tion defines implicitly the virial mass, M v i r , as that enclosed within r v - lT , 
and the virial velocity, V v i T , as the circular velocity measured at r v ; r . We 
compute A v j r (z) using A v i r (z) = 18tt 2 + 82/(2) - 39/(z) 2 , where 
Hz) = fno(l + z) 3 /(n o(l + Z ) 3 +^ A ))1-landn = fWl + ^bar 
('Bryan & Norman 1998), and is ~ 100 at z = 
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Figure 1. Star particles in one of our simulations, shown at z = 0. Particles are colored according to the age of the star; blue means a star is younger than 
~ 1 Gyr, red that it is older than ~ 10 Gyr. The large box is 2 r v j r (632 kpc) on a side and c entered on the primary galaxy. More than 85% of all stars are in 
the inner regions of the primary, within about ~ 20 kpc from the center (for more details see lAbadi et al.ll2006l) . surround the satellites "associated" with the 
primary galaxy; i.e., satellites that have been within r v j r in the past. Note that a few "associated" satellites lie well beyond the virial boundary of the system. 
Two of these satellites are highlighted for analysis in Figures|2]and|6] 



ISteinmetz & Navarrol d2002l) . The "primary" galaxies in these 
simulations have been analyzed in detail in several recent pa- 
pers, which the int er ested reader may wish to consult for detail s 
i lAbadi et alj2003allbtlMeza et al.l2003ll2005l ; lNavarro et alj2004h . 
We give a brief outline below for completeness. 

Each simulation follows the evolution of a small region of 
the universe chosen so as to encompass the mass of an L» galaxy 
system. This region is chosen from a large periodic box and res- 
imulated at higher resolution preserving the tidal fields from the 
whole box. The simulation includes the gravitational effects of dark 
matter, gas and stars, and follows the hydrodynamical evolution of 
the gaseous com ponent using the Smooth Particle Hydrodynamics 
(SPH) technique (Steinmetz 1996). We adopt the following cosmo- 
logical parameters for the ACDM scenario: Ho — 65 km/s/Mpc, 
a 8 = 0.9, fi A = 0.7, f2 C DM = 0.255, « b = 0.045, with no tilt in 
the primordial power spectrum. 

All re-simulations start at redshift zinit = 50, have force res- 
olution of order 1 kpc, and the mass resolution is chosen so that 
each galaxy is represented on average, at z = 0, with ~ 50, 000 
dark matter/gas particles. Gas is turned into stars a t rate s consistent 
with the empirical Schmidt-like law of lKennicuttl dl998h . Because 
of this, star formation proceeds efficiently only in high-density re- 
gions at the center of dark halos, and the stellar components of pri- 
mary and satellite galaxies are strongly segregated spatially from 
the dark matter. 

Each re-simulation follows a single ~ L* galaxy in detail, 
and resolves as well a number of smaller, self-bound systems of 
stars, gas, and dark matter we shall call generically "satellites". We 
shall hereafter refer to the main galaxy indistinctly as "primary" or 



"host". The resolved satellites span a range of luminosities, down 
to about six or seven magnitudes fainter than the primary. Each 
primary has on average ~ 10 satellites within the virial radius. 

Figure[T]illustrates the z — spatial configuration of star par- 
ticles in one of the simulations of our series. Only star particles are 
shown here, and are colored according to their age: stars younger 
than ~ 1 Gyr are shown in blue; those older than ~ 10 Gyr in red. 
The large box is centered on the primary and is 2 r v j r (632 kpc) on 
a side. The "primary" is situated at the center of the large box and 
contains most of the stars. Indeed, although not immediately ap- 
parent in this rendition, more than 85% of all stars are within ~ 20 
kpc from the center. Outside that radius most of the stars are old 
and belong to the stellar halo, except for a plume of younger stars 
stripped from a satellite that has recently merged with the primary. 
Satellites "associated" with the primary (see § l3,ll for a definition) 
are indicated with small boxes. Note that a few of them lie well 
beyond the virial radius of the primary. 

A preliminary analysis of the properties of the simulated satel- 
lite population and its relation to the stellar halo and the primary 
galaxy has been presented in Abadi, Navarro & Steinmetz (2006) 
and Sales et al (2007, submitted), where the interested reader may 
find further details. 
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Figure 2. Distance to the primary as a function of time for four satellites se- 
lected in one of our simulations. The four satellites are accreted into the pri- 
mary in two pairs of unequal mass. The heavier satellite of the pair, shown 
by solid lines, follows a "conventional" orbit: after turning around from the 
universal expansion, it accretes into the primary on a fairly eccentric orbit 
which becomes progressively more bound by the effects of dynamical fric- 
tion. Note that, once accreted, these satellites on "conventional" orbits do 
not leave the virial radius of the primary, which is shown by a dotted line. 
The light member of the pair, on the other hand, is ejected from the sys- 
tem as a result of a three-body interaction between the pair and the primary 
during first approach. One of the ejected satellites shown here is the "escap- 
ing" satellite identified in Figure[5] the other is the most distant "associated" 
satellite in that Figure. The latter is still moving toward apocenter at z = 0, 
which we estimate to be as far as ~ 3.5 r v ir- 



3 RESULTS AND DISCUSSION 
3.1 Satellites on conventional orbits 

The evolution of satellites in our simulations follows roughly the 
various stages anticipated by our discussion of the secondary infall 
model; after initially receding with the universal expansion, satel- 
lites turn around and are accreted into the primary. Satellites mas- 
sive enough to be well resolved in our simulations form stars ac- 
tively before accretion and, by the time they cross the virial radius 
of the primary, much of their baryonic component is in a tightly 
bound collection of stars at the center of their own dark matter ha- 
los. 

The stellar component of a satellite is thus quite resilient to 
the effect of tides and can survive as a self-bound entity for sev- 
eral orbits. This is illustrated by the solid lines in Figure [2] which 
show, for one of our simulations, the evolution of the distance to the 
primary of two satellites that turn around and are accreted into the 
primary at different times. As expected from the secondary infall 
model, satellites that are initially farther away turn around later; do 
so from larger radii; and are on more energetic orbits. After accre- 
tion (defined as the time when a satellite crosses the virial radius 
of the primary), their orbital energy and eccentricity are eroded by 
dynamical friction, and these two satellites do not leave the virial 
radius of the primary, shown by the dotted line in Figure [2] De- 
pending on their mass and orbital parameters, some of these satel- 



lites merge with the primary shortly after accretion, while others 
survive as self-bound entities until z = 0. For short, we shall refer 
to satellites that, by z — 0, have crossed the virial radius boundary 
at least once as satellites "associated" with the primary. 

The ensemble of surviving satellites at z — have kinematics 
consistent with the evolution described above. This is illustrated in 
Figure [3] where we show the radial velocities of all satellites as a 
function of their distance to the primary, scaled to virial units. Note 
that the majority of "associated" satellites (shown as circles in this 
figure) are confined within r v i r , and that their velocity distribution 
is reasonably symmetric and consistent with a Gaussian (Sales et 
al 2007). The most recently accreted satellites tend to have higher- 
than-average speed at all radii, as shown by the "crossed" circles, 
which identify all satellites accreted within the last 3 Gyr. 

Crosses (without circles) in this figure correspond to satellites 
that have not yet been accreted into the primary. These show a clear 
infall pattern outside r v i r , where the mean infall velocity decreases 
with radius and approaches zero at the current turnaround radius, 
located at about 3r v i r - All of these properties agree well with the 
expectations of the secondary infall model discussed above. 



3.2 Three-body interactions and satellites on unorthodox 
orbits 

Closer examination, however, shows a few surprises. To begin with, 
a number of "associated" satellites are found outside r v i r . As re- 



ported in_r£e\rous_work (see 1 _e ; g_ ; , Balogh et alj|2000l ; iMoore et al.l 
|2004| ; lGilletal.ll200l ; iDiemand et all 120071) . these are a minor- 
ity (~ 15% in our simulation series), and have been traditionally 
linked to departures from spherical symmetry during the accretion 
process. Indeed, anisotropies in the mass distribution during expan- 
sion and recollapse may endow some objects with a slight excess 
acceleration or, at times, may push satellites onto rather tangential 
orbits that "miss" the inner regions of the primary, where satellites 
are typically decelerated into orbits confined within the virial ra- 
dius. 

These effects may account for some of the associated satellites 
found outside r v i r at z = 0, but cannot explain why ~ 33% of all 
associated satellites are today on orbits whose apocenters exceed 
their turnaround radius. This is illustrated in Figure [4] where we 
show a histogram of the ratio between apocentric radius (measured 
at z — 0; r apo ) and turnaround radius (r ta ). The histogram high- 
lights the presence of two distinct populations: satellites on "con- 
ventional" orbits with r apo /r t a < 1, and satellites on orbital paths 
that lead them well beyond their original turnaround radius. 

Intriguingly, a small but significant fraction (~ 6%) of 
satellites have extremely large apocentric radius, exceeding their 
turnaround radius by 50% or more. These systems have clearly 
been affected by some mechanism that propelled them onto orbits 
substantially more energetic than the ones they had followed un- 
til turnaround. This mechanism seems to operate preferentially on 
low-mass satellites, as shown by the dashed histogram in Figure[4] 
which corresponds to satellites with stellar masses less than ~ 3% 
that of the primary. 

We highlight some of these objects in Figuref3] using "filled" 
circles to denote "associated" satellites whose apocenters at z = 
exceed their turnaround radii by at least 25%. Two such objects are 
worth noting in this figure: one of them is the farthest "associated" 
satellite, found at more than ~ 2.5 r v i r from the primary; the sec- 
ond is an outward-moving satellite just outside the virial radius but 
with radial velocity approaching ~ 2 VW- The latter, in particular, 
is an extraordinary object, since its radial velocity alone exceeds 
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Figure 3. Radial velocity of satellites versus distance to the primary. Ve- 
locities are scaled to the virial velocity of the system, distances to the virial 
radius. Circles denote "associated" satellites; i.e., those that have been in- 
side the virial radius of the primary at some earlier time. Crosses indicate 
satellites that are on their first approach, and have never been inside r v \ T . 
Filled circles indicate associated satellites whose apocentric radii exceed 
their turnaround radius by at least 25%, indicating that their orbital ener- 
gies have been substantially altered during their evolution. "Crossed" cir- 
cles correspond to associated satellites that have entered r v j r during the 
last 3 Gyrs. The curves delineating the top and bottom boundaries of the 
distribution show the escape velocity of an NFW halo with concentration 
c = 10 and c = 20, respectively. Note that there is one satellite "escap- 
ing" the system with positive radial velocity. Solid lines show the trajec- 
tories in the r — V r plane of the two "ejected" satellites shown in figure 
[2] Filled squares correspon d to th e fourteen brightest Milky Way satellites, 
taken from van den Bergl] j 19991) (complemented with NED data for the 
Phoenix, Tucana and NGC6822), and plotted assuming that V v ^ w - 109 
km/s and = 237 kpc (see Sales et al 2007). Arrows indicate how the 

positions of MW satellites in this plot would be altered if our estimate of 
V^. w (and, consequently, rM^) is allowed to vary by ±20%. 

the nominal escape velocit}0 at that radius. This satellite is on a 
trajectory which, for all practical purposes, will remove it from the 
vicinity of the primary and leave it wandering through intergalactic 
space. 

The origin of these unusual objects becomes clear when in- 
specting Figure [2] The two satellites in question are shown with 
dashed lines in this figure; each is a member of a bound pair of 
satellites (the other member of the pair is shown with solid lines of 
the same color). During first pericentric approach, the pair is dis- 
rupted by the tidal field of the primary and, while one member of 

2 The notion of binding energy and escape velocity is ill-defined in cosmol- 
ogy; note, for example, that the whole universe may be considered formally 
bound to any positive overdensity in an otherwise unperturbed Eistein-de 
Sitter universe. We use he re the nominal escape velocity of an NFW model 
iNavarro et aljl996t[l997l) to guide the interpretation. This profile fits rea- 
sonably well the mass distribution of the primaries inside the virial radius, 
and has a finite escape velocity despite its infinite mass. Certainly satellites 
with velocities exceeding the NFW escape velocity are likely to move far 
enough from the primary to be considered true escapers. 



the pair remains bound and follows the kind of "conventional" orbit 
described in § 13.11 the other one is ejected from the system on an 
extreme orbit. The trajectories of these two "ejected" satellites in 
the r-Vr plane are shown by the wiggly lines in Figure[3] 

These three-body interactions typically involve the first peri- 
centric approach of a bound pair of accreted satellites and tend to 
eject the lighter member of the pair: in the example of Figure [2] 
the "ejected" member makes up, respectively, only 3% and 6% of 
the total mass of the pair at the time of accretion. Other interac- 
tion configurations leading to ejection are possible, such as an un- 
related satellite that approaches the system during the late stages 
of a merger event, but they are rare, at least in our simulation se- 
ries. We emphasize that not all satellites that have gained energy 
during accretion leave the system; most are just put on orbits of 
unusually large apocenter but remain bound to the primary. This is 
shown by the filled circles in Figure [3j many affected satellites are 
today completing their second or, for some, third orbit around the 
primary. 

The ejection mechanism is perhaps best appreciated by in- 
specting the orbital paths of the satellite pairs. These are shown in 
Figure [6] where the top (bottom) panels correspond to the satellite 
pair accreted later (earlier) into the primary in Figure [2] Note that 
in both cases, as the pair approaches pericenter, the lighter member 
(dashed lines) is also in the process of approaching the pericenter 
of its own orbit around the heavier member of the pair. This co- 
incidence in orbital phase combines the gravitational attraction of 
the two more massive members of the trio of galaxies, leading to 
a substantial gain in orbital energy by the lightest satellite, effec- 
tively ejecting it from the system on an approximately radial orbit. 
The heavier member of the infalling pair, on the other hand, decays 
onto a much more tightly bound orbit. 

Figure [6] also illustrates the complexity of orbital configura- 
tions that are possible during these three-body interactions. Al- 
though the pair depicted in the top panels approaches the primary 
as a cohesive unit, at pericenter each satellite circles about the pri- 
mary in opposite directions: in the y-z projection the heavier mem- 
ber circles the primary clockwise whereas the ejected companion 
goes about it counterclockwise. After pericenter, not only do the 
orbits of each satellite have different period and energy, but they 
differ even in the sign of their orbital angular momentum. In this 
case it would clearly be very difficult to link the two satellites to 
a previously bound pair on the basis of observations of their orbits 
after pericenter. 

Although not all ejections are as complex as the one illustrated 
in the top panels of Figure [6] it should be clear from this figure 
that reconstructing the orbits of satellites that have been through 
pericenter is extremely difficult, both for satellites that are ejected 
as well as for those that remain bound. For example, the massive 
member of the late-accreting pair in Figure [2] sees its apocenter 
reduced by more than a factor of ~ 5 from its turnaround value 
in a single pericentric passage. Such dramatic variations in orbital 
energy are difficult to reproduce with simple analytic treatments in- 
spired on Chandrasekhar's dynamical friction formula (Penarrubia 
2007, private communication). 



4 APPLICATION TO THE LOCAL GROUP 

We may apply these results to the interpretation of kinematical out- 
liers within the satellite population around the Milky Way (MW) 
and M31, the giant spirals in the Local Group. Although part of the 
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Figure 4. Distribution of the ratio between the apocentric radius of satellites 
(measured at z = 0) and their turnaround radius, defined as the maximum 
distance to the primary before accretion. Note the presence of two groups. 
Satellites on "conventional" orbits have r apo /rt a < 1> the rest have been 
catapulted into high-energy orbits by three-body interactions during first 
approach. The satellite marked with a rightward arrow is the "escaping" 
satellite identified by a dot-centered circle in Figure[3] this system has nom- 
inally infinite r apo . The dashed histogram highlights the population of low- 
mass satellites; i.e., those with stellar masses at accretion time not exceeding 
2.6% of the primary's final M s tr- The satellite marked with an arrow is a 
formal "escaper" for which r apo cannot be computed. 



discussion that follows is slightly speculative due to lack of suit- 
able data on the three-dimensional orbits of nearby satellites, we 
feel that it is important to highlight the role that the concomitant 
accretion of multiple satellites may have played in shaping the dy- 
namics of the dwarf members in the Local Group. 



4.1 Milky Way satellites 

The filled squares in Figure[3]show the galactocentric radial veloc- 
ity of thirteen bright satellites around the Milky Way and compare 
them with the simulated satellite population. This comparison re- 
quires a choice for the virial radius and virial velocity of the Milky 
Way, which are observationally poorly constrained. 

We follow here the approach of Sales et al (2007), and use the 
kinematics of the satellite population itself to set the parameters of 
the Milky Way halo. These authors find that simulated satellites 
are only mildly biased in velocity relative to the dominant dark 
matter component: a r ~ 0.9(±0.2)VW, where cr r is the radial 
velocity dispersion of the satellite population within r V i r . Using 
this, we find = 109 ± 22 km/s and r^ r w = 237 ± 50 kpc 

from the observed radial velocity dispersion of ~ 98 km/s. This 
corresponds to Mjf r w = 7 x 10 11 Mp>, in reasonable a greement 



with the 1-2 x 10 12 M^ estimate of 



Klyp in et alj d2002l) and with 
the recent findings of lSmith et al.1 d2006h based on estimates of the 
escape velocity in the solar neighbourhood. 

Since Leo I dwarf has the largest radial velocity of the Milky 
Way satellites, we have recomputed the radial velocity dispersion 



excluding it from the sample. We have found that ay drops from 
98 to 82 km/s when Leo I is not taken into account changing our 
estimation of V V I ^ W from 109 to 91 km/s, still within the errors 
of the value previously found. Given the recent rapid growth in 
the number of known Milky Way satellite one would suspect that 
the velocity dispersion will significantly increase if more Leo I- 
like satellites are detected. However, we notice that given their high 
velocities they are not expected to remain inside the virial radius for 
a long time period hence not contributing to the ay computation. 

Figure [3] shows that, considering V^. w = 109 km/s, the ve- 
locities and positions of all MW satellites are reasonably consis- 
tent with the simulated satellite population, with the possible ex- 
ception of Leo I, which is located near the virial radius and is 
moving outward with a velocity clearly exceeding V v i r . Indeed, for 



' VI 



109 km/s, Leo I lies right on the escape velocity curve 



of an NFW profile with concentration parameter similar to those 
measured in the simulations. This is clearly a kinematical outlier 
reminiscent of the satellite expelled by three-body interactions dis- 
cussed in the previous subsection and identified by a dot-centered 
circle in Figure[3] This is the only "associated" satellite in our simu- 
lations with radial velocity exceeding Kdr and located outside r v ir- 
Could Leo I be a satellite that has been propelled into a highly- 
energetic orbit through a three-body interaction? If so, there are a 
number of generic predictions that might be possible to verify ob- 
servationally. One is that its orbit must be now basically radial in 
the rest frame of the Galaxy, although it might be some time before 
proper motion studies are able to falsify this prediction. A second 
possibility is to try and identify the second member of the pair to 
which it belonged. An outward moving satellite on a radial orbit 
takes only ~ 2-3 Gyr to reach r v i r with escape velocity. Coinciden- 
tally, this is about the time that the Magellanic Clouds pair were last 
at pericenter, acc ording to the traditional orbital evolution of the 
Cloud s (see, e.g., I Gar diner & Noguchll ll996l : Ivan der Mar el et al.1 
|2002|) . 

Could Leo I have been a Magellanic Cloud satellite ejected 
from the Galaxy a few Gyrs ago? Since most satellites that are 
ejected do so during first pericentric approach, this would im- 
ply that the Clouds were accreted only recently into the Galaxy, 
so that they reached their first pericentric approach just a few 
Gyr ago. This is certainly in the spirit of the re-analysis of the 
orbit of the Clouds presented recently by iBesla et aU |2007| and 
based on new proper motion measurements recently reported by 
iKallivavalil etafl hoOfj) . In this regard, the orbit of the Clouds 
might resemble the orbit of the companion of the "escaping" satel- 
lite located next to Leo I in Figure [3] The companion is fairly 
massive and, despite a turnaround radius of almost ~ 600 kpc 
and a rather late accretion time (t acc = 10.5 Gyr, see Figure^, 
it is left after pericenter on a tightl y bound, short-period orbit 
resembling that of the Cl ouds today ( Gar diner & Noguchi|[T99r3 : 
van der Marel et al. 2002). To compound the resemblance, this 
satellite has, at accretion time, a total luminosity of order ~ 10% 
of that of the primary, again on a par with the Clouds. 

We also note that an ejected satellite is likely to have picked 
up its extra orbital energy through a rather close pericentric passage 
and that this may have led to substantial tidal damage. This, indeed, 
has been argued recently bv lSohn et alj2006l on the basis of asym- 
me tries in the spatia l and velocity distribution of Leo I giants (but 
see lKoch et aU2007l for a radically different interpretation). 

On a final note, one should not forget to mention another (less 
exciting!) explanation for Leo I: that our estimate of is a sub- 
stantial underestimate of the true virial velocity of the Milky Way. 
The arrows in Figure[3]indicate how the position of the MW satel- 
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Figure 5. As Figure [3] but for line-of-sight velocities and projected dis- 
tances. Three random orthogonal projections have been chosen for each 
simulated satellite system. Signs for Vi os have been chosen so that it is 
positive if the satellite is receding away from the primary in projection, 
negative otherwise. The "escaping" satellite from Figure [5] is shown by a 
starred symbol. Filled square s correspond t o the M3 1 satellites ta ken from 
iMcConnachie & Irwinl200d plus And XIV <Maiewski et al.l2007t) and And 
XII (Chapman et al 2007, submitted) and assuming that V^ 31 ~ 138 km/s 
and r^ 31 = 300 kpc. Arrows indicate how the positions of M31 satellites 
in this plot would be altered if our estimate of V^ 31 (and, consequently, 
r^ 31 ) is allowed to vary by 20%. 

lites in this plane would change if our estimate of V V I ^ W is varied by 
±20%. Increasing V^ W by ~ 20% or more would make Leo I's 
kinematics less extreme, and closer to what would be expected for a 
high-speed satellite completing its first orbit. This rather more pro- 
saic scenario certainly cannot be discounted on the basis of avail- 
able data (see, e.g., Zaritsky et al 1989, Kochanek 1996, Wilkinson 
& Evans 1999) 



4.2 M31 satellites 

A similar analysis may be applied to M31 by using the projected 
distances and line-of-sight velocities of simulated satellites, shown 
in Figure[5] Three orthogonal projections of the simulated satellites 
are overlapped in this figure, with symbols as defined in Figure[3] 
Following the same approach as in § 14, 11 we use the fact that the 
line-of-sight satellite velocity dispersion is <ti os ~ 0.8(±0.2) V v i r 
in our simulations to guide our choice of virial velocity and radius 
for M31; V™ 31 = 138 ± 35 km/s and r^ 31 = 300 ± 76 kpc. 
(We obtain a\ os = 111 km/s for all 17 satellites within 300 kpc of 
M31.) This compares favo urably with the Vffi 31 ~ 120 km/s esti- 
mate recently obtained bv lSeigar et al. 2006 under rather different 
assumptions. 

With this choice, we show the 19 satellites around M31 
compiled by McConnachie & Irwin (2006), plus two recently- 
discovered satellites for which positions and radial velocities have 




-50 50 -50 50 -50 50 
x (kpc) y (kpc) z (kpc) 




x (kpc) y (kpc) z (kpc) 

Figure 6. Orbital paths for both pair of satellites shown in Figure [2] Upper 
(bottom) panel shows the pair that accretes later (earlier) in that figure and 
shows the orbits in the rest frame of the primary. The coordinate system is 
chosen so that the angular momentum of the primary is aligned with the 
z axis. A solid curve tracks the path of the heavier satellite; a dashed line 
follows the satellite that is propelled into a highly energetic orbit after. 



become available (And XII, Chapman et al 2007, and And XIV, 
Majewski et al. 2007). As in Figure[3] arrows indicate how the po- 
sition of M3 1 satellites would change in this figure if V^ 31 were 
allowed to vary by ±20%. We notice that the exclusion of And XII 
and And XIV (the highest velocity satellites within 300 kpc from 
Andromeda) in the V^ 31 estimation gives ~ 100 km/s, consistent 
with the V^ 31 = 138 ± 35 km/s previously found considering all 
satellites. Projected distances are as if viewed from infinity along 
the direction joining the Milky Way with M31 and that the sign of 
the line-of-sight velocity in Figure|5]is chosen to be positive if the 
satellite is receding from the primary (in projection) and negative 
otherwise. 

There are a few possible outliers in the distribution of M31 
satellite velocities: And XIV ( Majewski et al 2007), the Pega- 
sus dwarf irregular (UGC 12613, iGallagher et alj[l998l) . And XII 
(Chap man et al 2007), and UGCA 092 (labelled U092 in Fig- 
ure El M cConnachie & Irwinl kooe). And XIV and PegDIG seem 
likely candidates for the three-body "ejection" mechanism dis- 
cussed above: they have large velocities for their position, and, 
most importantly, they are receding from M31; a requirement for 
an escaping satellite. Note, for example, that And XIV lies very 
close to the "escaping" satellite (dot-centered symbol in Figure [5} 
paired to Leo I in the previous subsection. Escapers should move 
radially away from the primary, and they would be much harder 
to detect in projection as extreme velocity objects, unless they are 
moving preferentially along the line of sight. It is difficult to make 
this statement more conclusive without further knowledge of the 
orbital paths of these satellites. Here, we just note, in agreement 
with Majewski et al (2007), that whether And XIV and PegDIG are 
dynamical "rogues" depends not only on the (unknown) transverse 
velocity of these galaxies, but also on what is assumed for M31's 
virial velocity. With our assumed V^ 31 = 138 km/s, neither And 
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XIV nor PegDIG look completely out of place in Figure [5j had 
we assumed the lower value of 120 km/s advocated by Seigar et 
al (2006) And XIV would be almost on the NFW escape velocity 
curve, and would certainly be a true outlier. 

High- velocity satellites approaching M31 in projection are 
unlikely to be escapers, but rather satellites on their first approach. 
This interpretation is probably the most appropriate for And XII 
and UGCA 092. As discussed by Chapman et al (2007), And XII 
is almost certainly farther than M31 but is approaching us at much 
higher speed (~ 281 km/s faster) than M31. This implies that And 
XII is actually getting closer in projection to M3 1 (hence the nega- 
tive sign assigned to its V\ oa in Figure[5j, making the interpretation 
of this satellite as an escaping system rather unlikely. 

Note, again, that although And XII (and UGCA 092) are just 
outside the loci delineated by simulated satellites in Figure [5] re- 
vising our assumption for V^ 31 upward by 20% or more would 
render the velocity of this satellite rather less extreme, and would 
make it consistent with that of a satellite on its first approach to 
M3 1 . As was the case for Leo I, this more prosaic interpretation of 
the data is certainly consistent with available data. 



5 SUMMARY AND CONCLUSIONS 

We examine the orbits of satellite galaxies in a series of 
Nbody/gasdynamical simulations of the formation of L» galax- 
ies in a ACDM universe. Most satellites follow orbits roughly in 
accord with the expectations of secondary infall-motivated mod- 
els. Satellites initially follow the universal expansion before being 
decelerated by the gravitational pull of the main galaxy, turning 
around and accreting onto the main galaxy. Their apocentric radii 
decrease steadily afterwards as a result of the mixing associated 
with the virialization process as well as of dynamical friction. At 
z = most satellites associated with the primary are found within 
its virial radius, and show little spatial or kinematic bias relative to 
the dark matter component (see also Sales et al 2007). 

A number of satellites, however, are on rather unorthodox or- 
bits, with present apocentric radii exceeding their turnaround radii, 
at times by a large factor. The apocenters of these satellites are 
typically beyond the virial radius of the primary; one satellite is 
formally "unbound", whereas another is on an extreme orbit and is 
found today more than 2.5r v j r away, or > 600 Mpc when scaling 
this result to the Milky Way. 

These satellites owe their extreme orbits to three-body inter- 
actions during first approach: they are typically the lighter member 
of a pair of satellites that is disrupted during their first encounter 
with the primary. This process has affected a significant fraction of 
satellites: a full one-third of the simulated satellite population iden- 
tified at z — have apocentric radii exceeding their turnaround 
radii. These satellites make up the majority (63%) of systems on 
orbits that venture outside the virial radius. 

We speculate that some of the kinematical outliers in the Local 
Group may have been affected by such process. In particular, Leo 
I might have been ejected 2-3 Gyr ago, perhaps as a result of in- 
teractions with the Milky Way and the Magellanic Clouds. Other 
satellites on extreme orbits in the Local Group may have origi- 
nated from such mech anism. Cetus dLewis etaU 2007) and Tucana 
dOosterloo et al J f 1996) — two dwarf spheroidals in the periphery 
of the Local Group — may owe their odd location (most dSphs are 
found much closer to either M31 or the Galaxy) to such ejection 
mechanism. 

If this is correct, the most obvious culprits for such ejection 



events are likely to be the largest satellites in the Local Group (M33 
and the LMC/SMC), implying that their possible role in shaping 
the kinematics of the Local Group satellite population should be 
recognized and properly assessed. In this regard, the presence of 
kinematical oddities in the population of M31 satellites, such as 
the fact that the majority of them lie on "one side" of M31 and 
seem to be receding away from it (McConnachie & Irwin 2006), 
suggest the possibility that at least some of the satellites normally 
associated with M3 1 might have actually been brought into the Lo- 
cal Group fairly recently by M33. Note, for example, that two of 
the dynamical outliers singled out in our discussion above (And 
XII and And XIV) are close to each other in projection; have rather 
similar line-of-sight velocities (in the heliocentric frame And XII 
is approaching us at 556 km/s, And XIV at 478 km/s); and belong 
to a small subsystem of satellites located fairly close to M33. 

The same mechanism might explain why the spatial distribu- 
tion of at least some satellites, both around M3 1 and the Milky Way, 
seem to align themselves on a "planar" configur ation jMa iewski 
ll994l ; lLibeskind et alj200H : lKoch & Grebell2006h . as this may just 
reflect the orbital accretion plane of a multi ple system of satel- 
lites accreted sim ultaneously in the recent past ( Krou paet alj2005l : 
iMetz et all2007r) . 

From the point of view of hierarchical galaxy formation mod- 
els, it would be rather unlikely for a galaxy as bright as M33 to 
form in isolation and to accrete as a single entity onto M31. There- 
fore, the task of finding out which satellites (rather than whether) 
have been contributed by the lesser members of the Local Group, 
as well as what dynamical consequences this may entail, should 
be undertaken seriously, especially now, as new surveys begin to 
bridge our incomplete knowledge of the faint satellites orbiting our 
own backyard. 
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